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Vector  stochastic  variational  principles  are  derived  for  the  statistics 
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20.  the  scatterers.  The  nonstochastic  nature  of  the  Incident  field  allows 
the  statistical  moments  of  T  and  of  the  differential  scattering  cross 
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(4ir)-n<NinN2n/Dn)»  and  should  allow  practical  application  of  variational 
techniques  to  random  scattering  problems. 
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COOPERATIVE  EFFECTS  IN  LIGHT  SCATTERING  FROM  AEROSOL 
CLUSTERS  AND  SURFACE  ASSEMBLIES 


1 .  INTRODUCTION 

1.1  Objectives 

The  long-range  objective  of  this  program  Is  to  acquire  a  better 
understanding  of  cooperative  effects,  l.e..  Interference  and  multiple 
scattering,  and  polarization  In  the  scattering  and  absorption  of  electro¬ 
magnetic  radiation  by  random  surfaces  and  random  collections  of  particles. 
This  objective  Is  to  be  achieved  by  developing  stochastic  variational 
principles  for  calculating  wave  scattering  and  using  specialized- trial 
cases  to  test  the  efficacy  of  these  principles. 

1.2  Progress 

During  the  short  span  of  the  reporting  period,  we  have  completed 
the  development  of  vector  stochastic  variational  principles  for  electro¬ 
magnetic  plane  wave  scattering  from  perfectly  conducting  objects  and  rough 
surfaces  and  from  conducting  dielectrics  which  are  generally  inhomogeneous 
and  anisotropic.  The  parameters  describing  the  scatterers  may  be  deter¬ 
ministic  or  random  with  completely  arbitrary  statistics.  Specifically, 
the  stochastic  vector  principle  for  conducting  dielectrics  was  extended 
to  Include  anisotropic  materials  with  no  symmetry  restrictions  on  the 
conductivity  and  permittivity  tensors.  In  addition,  a  stochastic  vari¬ 
ational  principle  was  derived  for  polarized  plane  wave  scattering  from 
a  randomly  rough  planar  surface  which  Is  perfectly  conducting.  This 
formulation  incorporates  a  half-space  dyadic  Green  function.  Finally,  a 
variational  principle  for  closed,  perfectly  conducting  scatterers  was 
recast  Into  stochastic  form.  These  variational  principles  are  for  the 
statistical  moments  and  probability  density  functions  of  T  and  |T|2,  where 

T  =  es*T  Is  the  component  of  the  normalized  far-field  vector  scattering 

amplitude  T  along  with  an  arbitrary  polarization  direction  e%.  A  detailed 

description  of  this  work  Is  given  In  Sections  2  through  6  which  follow. 

These  sections  are  an  exact  duplication  of  a  report  by  Jerry  A.  Krill 
and  Robert  H.  Andreo  of  the  Applied  Physics  Laboratory,  The  Johns 
Hopkins  University. 

1.3  Future  Work 

Our  work  In  the  near  term  will  Involve  the  development  of  a  completely 
general  stochastic  vector  variational  principle  for  scatterers  with  in¬ 
homogeneous  and  anisotropic  magnetic,  as  well  as  dielectric  and  conductive, 
properties.  Subsequent  efforts  will  then  be  directed  towards  variational 
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calculations  for  specific  scattering  problems.  This  work  will  be  supported 
by  the  Army  Research  Office. 

2.  VECTOR  STOCHASTIC  VARIATIONAL  PRINCIPLES  FOR  ELECTROMAGNETIC  WAVE 

SCATTERING 

This  Is  a  report  by  J.  A.  Krill  and  R.  H.  Andreo  of  The  Johns  Hopkins 
University,  Applied  Physics  Laboratory,  Johns  Hopkins  Road,  Laurel,  Mary¬ 
land  20810. 

2.1  History 

During  the  past  three  decades  numerous  approximation  techniques  have 
been  applied  to  model  electromagnetic  wave  scattering  from  surfaces  and 
media.  ~3  Variational  methods4-6  have  been  of  value  in  studies  of  deter¬ 
ministic  scattering  and  propagation  and  are  recognized  as  having  several 
important  advantages.4-1®  The  primary  one  is  the  cancellation  of  first 
order  errors  arising  from  the  initial  estimate  of  the  (generally  unknown) 
field  on  or  within  the  scatterer.  This  cancellation  will  lead  to  accurate 
results  if  the  estimate  is  already  a  good  approximation.  Thus,  it  is  for¬ 
tunate  that  the  variational  principle  for  scattering  contains  a  built-in 
indicator  of  accuracy  and  region  of  validity.7  Also,  in  certain  cases  an 
error  analysis  may  be  performed  which  provides  a  useful  error  bound  for 
the  chosen  trial  estimate.11  Until  recently,7  however,  difficulties  in¬ 
herent  in  calculating  the  statistics  of  the  scattered  field  using  the 
variational  principle  have  precluded  more  than  limited  application  to 
stochastic  problems.7*8*11  The  difficulties  occur  in  evaluating  the 
statistical  moments  of  the  ratio  of  complex  integrals,  a  characteristic 
of  variational  scattering  calculations. 

In  1977,  Hart  and  Farrell7  recast  a  scalar  variational  principle 
into  a  form  which  does  not  possess  this  inherent  complexity,  thus  allowing 
tractable  variational  evaluation  of  the  mean  of  the  scattering  amplitude 
T  and  of  its  absolute  square  |T|2.  Requiring  no  assumptions  concerning 
the  statistical  nature  of  the  scatterer,  they  proved  that  the  average  values 
of  the  individual  integrals  comprising  a  deterministic  variational  expres¬ 
sion  (e.g.,  for  T)  may  themselves  be  combined  to  form  a  variationally  in¬ 
variant  expression.  Because  evaluation  of  the  mean  of  each  integral  is 
inherently  more  tractable  than  evaluation  of  the  mean  of  the  integral  com¬ 
bination,  this  new  variational  approach  is  expected  to  allow  more  tractable 
calculation  of  the  statistics  of  the  scattered  field.  It  has  therefore 
been  labeled  the  "stochastic"  variational  principle.12 

Gray,  Hart,  and  Farrell13*14  performed  the  first  application  of  this 
stochastic  principle  by  treating  a  perfectly-conducting  random  rough  sur¬ 
face  model.  This  surface  consists  of  parallel,  infinitely-long  hemicylin- 
ders  randomly  distributed  over  an  infinite  plane  but  constrained  to  be 
nonoverlapping.  A  plane  electromagnetic  wave  is  incident  normally  to  and 
polarized  parallel  to  the  hemicylinder  axes,  thus  reducing  this  problem  to 
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a  scalar,  two-dimensional  treatment.  In  the  limit  of  a  large  number  of 
Rayleigh  hemlcyllnders  they  obtained  the  Born  (l.e.,  the  first  order  per¬ 
turbation)  approximation15  for  the  mean  absolute-squared  scattering  ampli¬ 
tude  (|T|5)  In  closed,  analytic  form.  They  then  performed  a  variational 
Improvement  of  the  Born  approximation,  obtaining  a  closed-form  variational 
solution,  and  noted  a  discrepancy  between  these  results.  This  discrepancy 
was  studied  further  by  Krill  and  Farrell16  who  examined  a  special  case  of 
this  surface  In  which  only  two  Rayleigh  hemlcylinders  are  present.  By 
comparing  the  exact  solution  to  the  corresponding  Born  and  variational 
approximations  of  (|T|2>  for  this  surface,  they  demonstrated  that  the  var¬ 
iational  approximation  is  the  more  accurate  one.  This  greater  accuracy 
is  achieved  because  the  variational  result  accounts  for  multiple  scatter¬ 
ing,  which  Is  significant  in  this  case,  whereas  the  Born  approximation 
does  not. 

Thus  far,  the  stochastic  variational  expressions  and  their  applica¬ 
tions  have  been  appropriate  either  for  scalar  (e.g.,  acoustic)  fields  or 
for  vector  fields  with  special  scattering  configurations.13’16  Also,  they 
have  involved  only  homogeneous  boundary  conditions  (e.g.,  perfect  conduc¬ 
tors).  Vector  stochastic  variational  principles  will  be  derived  in  this 
paper  to  account  for  arbitrary  electromagnetic  field  polarization  and  to 
allow  treatment  of  conducting  and  dielectric  random  scatterers,  which  may 
be  inhomogeneous  and  anisotropic.  These  new  stochastic  variational  formu¬ 
lations  express  the  statistical  moments  of  the  components  of  the  vector 

scattering  amplitude  T  as  quotients  of  simpler  statistical  moments  in  a 
manner  analogous  to  the  scalar  development  of  Hart  and  Farrell.7 

2.2  Scope  of  This  Treatment 

As  in  the  scalar  case,  the  vector  stochastic  variational  expressions 
are  based  on  corresponding  deterministic  ones.  Hence,  in  sections  3  and  4 
we  present  deterministic  variational  principles  for  vector  wave  scattering 
from  arbitrary  conducting  dielectrics,  from  a  perfectly-conducting,  rough 
planar  surface,  and  from  arbitrary  perfect  conductors.  The  formulations 
for  the  two  configurations  involving  perfect  conductors  are  distinguished 
by  the  choice  of  dyadic  Green  function.  The  vector  variational  formulations 
for  scattering  from  arbitrary,  nonplanar  perfect  conductor  has  appeared 
previously10  and  is  included  here  for  completeness.  In  section  5  we  prove 

the  variational  invariance  of  the  expressions  for  the  arbitrary  n^  statis- 

tlcal  moments  of  the  vector  components  of  T,  and  of  their  absolute  squares, 
based  on  the  deterministic  formulations  of  the  previous  sections.  Finally 
section  6  discusses  the  results  and  indicates  some  potential  physical 
applications. 

3.  VARIATIONAL  PRINCIPLE  FOR  CONDUCTING,  DIELECTRIC  SCATTERERS 

In  this  section  we  shall  derive  in  some  detail  a  variational  expres¬ 
sion  for  the  scattering  of  a  plane  electromagnetic  wave  by  a  finitely  con- 
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ducting,  dielectric  volume,  based  on  the  vector  wave  equation  for  the 

electric  field  E.  It  Is  an  adaptation,  to  (vector)  electromagnetic  waves, 
of  the  scalar  variational  principle  originally  developed  by  Schwinger  for 
the  Schrodlnger  wave  equation  In  quantum  mechanics,15  and,  to  our  knowledge, 
has  not  appeared  previously.  (See,  however,  References  9  and  17-19  for 
related,  but  somewhat  different,  variational  principles.)  This  expression 
will  be  recast  Into  a  stochastic  variational  principle  in  section  5. 

Figure  1  Illustrates  the  scattering  problem.  A  plane  wave  with 

*>■ 

i  lc*f  *  r 

electric  field  E'(r)  =  e  is  Incident  upon  a  scatterer  (or  scat¬ 
tered)  with  free  space  permeability  uo  and  with  generally  inhomogeneous, 

anisotropic  (i.e.,  dyadic)  conductivity  c(r)  and  permittivity  ¥(r) .  The 
scattering  volume  Vo  is  suspended  in  free  space  (e0,  y0)  and  its  geometri¬ 
cal  characteristics,  such  as  shape  and  orientation,  are  arbitrary.  We 
wish  to  examine  a  vector  polarization  component  of  the  scattering  amplitude 

T  along  an  arbitrary  direction  e  ,  i.e.,  es*T,  at  an  arbitrary  position 

r  in  the  far-fleld  (r  -*■  «,  r  =  |r|). 

The  time-harmonic  electric  field  E  is  governed  by  the  vector  wave 
equation,10’20 

V*VxE(r)  -  kjE(r)  =  U(r)*E(r),  (1) 

for  which  we  have  defined 


U(r)  =  k2  1  ’  (2) 

where  k2  =  oj2poeo,  w  is  the  angular  frequency,  and  I  =  xx  +  yy  +  zz  is  the 
unit  dyadic.  The  total  field  E  can  be  expressed  as  the  sum  of  the  incident 

-*■  -*■•5  +5C  +.(■ 

and  scattered  contributions,  i.e.,  E  =  E  +  E  ,  and  E  satisfies  the 
usual  radiation  condition  at  infinity,20 


lim  r 

r  +  oo 


JvxEsc(r)  "  rxEsc(r^|  = 


for  any  localized  current  distribution  on  the  scatterer,  where  r  =  r/r. 

■f 

This  allows  the  asymptotic  form  (r  »  r')  of  E  to  be  written  as9*10 


ik.*r 


E(r)  >  Aai  e  1  +  AT(ks,ki)e1k°r/r 
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for  any  localized  scatterer.  This  equation  defines  the  (angular-dependent) 
scattering  amplitude  T(k  ,k.).  Me  require  T  to  be  perpendicular  to  the 

■4*  S  1 

scattered  wave  vector  k  since  the  scattered  wave  In  the  far  field  Is 
transverse.10*20 

4 

In  order  to  obtain  an  Integral  equation  for  E  based  on  (1),  we  apply 
Green's  Theorem  for  the  vector  wave  equation27  to  yield 


-4  -4  -44  -4 

E(r)  =  E  (r)  + 


|  G0(r,r')-p(r*  )-E(r*  rjdV* , 


where  the  volume  V  Is  taken  to  be  all  of  space,  bounded  only  by  the  surface 

at  Infinity,  on  which  the  surface  Integrals  involving  E  vanish  by  virtue 

of  the  radiation  condition.  The  work  of  Van  Blade!27  and  Yaghjlan28  showed 

that  (5)  Is  valid  even  In  the  "source"  region  U(r')  f  0,  provided  that  one 
chooses  the  dyadic  Green  function 

G„(r,r')  =  P.V.  G(r,r')  +  I6(r-r'),  (6) 

-  — ^  ^  ■  ■"  i  ^ 

where  G(r,r')  «  (I+W/k2)exp(ik0R)/4irR  is  the  familiar  free  space  dyadic 

Green  function,20  with  R  =  |r-r'|.  They  demonstrated  this  result  by 
applying  Green's  Theorem  while  carefully  excluding  the  volume  around  the 

singular  point  r'  =  r.  The  symbol  P.V.  In  (6)  indicates  that  the  principle 

value  Is  to  be  taken  of  an  Integral,  such  as  that  in  (5),  which  involves  G. 
The  second  term  In  (6)  arises  because  of  the  nature  of  the  singularity  in 

the  term  W(l/R)  at  R  =  0,  and  the  dyadic  L  depends  on  the  shape  of  the 
excluded  volume  used  In  computing  the  principal  value  integral.  Explicit 

expressions  for  L  for  the  more  common  excluded  volume  shapes  have  been 

tabulated  In  Reference  28.  Note  that  Go  reduces  to  G  when  r  is  outside 

the  source  region;  however,  the  more  general  form  G0  will  be  needed  In 
the  development  of  a  variational  principle. 

Applying  the  asymptotic  form  of  G0  to  (5)  and  observing  (4)  leads 
to  an  Integral  expression  for  T, 


Ni (ks,ki) 


(8) 


where  k  h  kor  and  =  I  -  rr  is  the  projection  operator  onto  the  plane 
s  -v  r 

transverse  to  k$. 

-► 

To  complete  our  derivation  of  the  variational  expression  for  T,  we 
first  premultiply  (5)  by  a  vector  [U(r)*E(r)]^r,  where  the  superscript 
Tr  denotes  matrix  transposition.  The  operator  U  and  the  field  E  are  as 
yet  undetermined  except  that,  analogous  to  U,  we  choose  U  to  be  independent 

-f 

of  E  and  E.  The  result  is  integrated  over  rcV  to  obtain 


A  =  D(ks,ki)/N2(ks,k1) 


(9) 


with 


Mks,k.) 


U(r)-E(r) 

Tr 

■ 

ik.*r 

A  _  1 

e.  e 

V 

- 

. 

(10) 


and 


D<ks*ki> 


V  L 


-.Tr 


U(r)*E(r) 


■E(r)dV 


■ 

p 

— —  ^ 

rTr 

—  ->■ 

dv 

dV' 

U(r)-E(r) 

•Go(r.r')* 

U(r')*E(r') 

V 

V 

_  _ 

- 

The  evaluation  of  the  double  integral  in  (11)  requires  the  dyadic  Green 
function  in  the  source  region,  and  the  definition  given  in  (6)  guarantees 
that  this  double  integral  is  well  defined.  We  combine  (7)  and  (9)  and 
premultiply  the  result  by  an  arbitrary,  constant  vector  e$  to  yield 

T(k$,k.)  =  ±  N1(ks,k.)N2(ks,ki)/D(ks,k.),  (12) 
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■4  ■+ 

where  Ni  =  £s*Nj,  and  where  T  h  es*T  is  the  (experimentally  measurable) 
component  of  the  vector  scattering  amplitude  along  e„.  We  note  that 

5  4- 

expressions  (7),  (9),  and  (12)  are  equalities  only  if  the  exact  field  E, 
l.e.,  the  solution  to  (5),  is  used  in  the  evaluation  of  Ni,  N2,  and  D. 

We  will  now  examine  what  happens  when  approximate  fields  (which  do 
not  in  general  satisfy  the  wave  equation  and  boundary  conditions)  are  used. 

It  will  be  shown  that,  for  specific  conditions  on  E,  expression  (12)  be¬ 
comes  variationally  Invariant  in  that  the  first-order  variation  of  the 

functional  T  about  the  exact  fields  E  and  E  vanishes, 9»21"2*  i.e.,  <5T  =  0. 
This  implies  that  first-order  errors  in  Ni,  N2,  and  D  which  arise  from 

*f  4  4  4 

replacing  the  exact  fields  E  and  E  in  (12)  by  approximate  fields  E  +  6E 

-4  -4 

and  E  +  6E  do  not  lead  to  first  order  errors  in  the  corresponding  approxi¬ 
mation  for  T.  We  observe  that,  upon  using  the  same  approximate  field 

-4  -4 

E  +  (5E  to  evaluate  Nx,  (7)  also  becomes  an  approximation  to  T,  but  with 
first-order  errors.*  To  determine  the  conditions  for  invariance,  the 

-4  -4  -4  -4 

approximate^fields  are  first  inserted  into  (12)  to  obtain  T[E  +  6E,  E  +  6E], 

*4  4 

where  T[E,_E]  denotes  the  exact  expression.  Expanding  the  difference 

TIE  +  6E,  E  +  6E]  -  T[E,  E]  in  terms  of  the  arbitrary  variations  6E  and 

6E,  retaining  only  first-order  terms,  and  dividing  the  result  by  the 
exact  expression  (12),  one  may  show  that  the  variation21  6T  is 

6T/T  =  6Ni/Nj  +  5N2/N2  -  6D/D,  (13) 


where  6Ni,  6N2,  and  6D  are  the  variations  of  Ni,  N2,  and  D,  respectively. 

-4T  v* 

The  condition  6T  =  0  and  the  arbitrariness  of  the  variation  6E  in 
(13)  identifies  the  Euler  equation20’22  for  E, 


->  ik,*r 

E(r)  =  Ae.  e  1  + 


dV'  G0(r,r')*|u(r,)»E(r,)J,  reV„  (14) 


and  the  arbitrariness  of  the  variation  6E  identifies  (after  taking  the 
transpose)  the  Euler  equation  for  E, 


♦Equations  (7)  and  (12)  lead  to  different  approximations  for  T  because  the 
ratio  (D/N2)  is  not  equal  to  the  amplitude  of  the  Incident  wave  when  approx¬ 
imate  fields  are  used. 
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dV  G0Tr(r,r')*U(r)-E(r),  r'eV0 

(15) 

where  the  superscript  -1  denotes  matrix  inversion.  Observe  that  the  Euler 

-►  -►  = 
equations  determine  E  and  E  only  within  the  volume  Vo  where  U  is  non-zero, 
i.e.,  the  volume  of  the  scatterer.  (Of  course,  V0  may  coincide  with  V, 
e.g.,  for  a  scatterer  with  "soft"  boundaries.)  Thus,  the  Euler  equation 

(14)  for  E  Is  identical  to  the  wave  equation  (5)  only  within  Vo.  This  is 
consistent  with  the  observation  that  the  scattering  amplitude  (7)  is  com¬ 
pletely  determined  by  the  field  within  the  scatterer. 

Conversely,  provided  solutions  for  E  and  E  exist,22’23  these  Euler 

equations  imply  the  vanishing  of  6T.  Thus,  under  those  conditions 

which  the  Euler  equations  may  be  satisfied,  (12)  is  variationally  invariant, 

“►  •+ 

in  that  first  order  variations  about  the  exact  solutions  E  and  E  are  zero. 

(See  References  22  and  23  for  further  discussion  of  these  conditions.) 

For  our  purposes  it  is  convenient  to  choose  the  undetermined  operator 

=  =  =  =rr  ■* 

U  to  equal  the  transpose  of  U  so  that  U  =  U  .  The  Euler  equation  for  E 
then  reduces  to 


U(r' ) *E(r ' ) 


N7  V. 


•Ik  *r' 


W*)  -ftVs 


-ik$-r' 


dV  G0(r',r) 


(16) 


for  r‘eV0,  where  we  have  used  the  symmetry  properties  of  G0.20  This  equa¬ 
tion  is  now  in  the  form  of  an  integral  wave  equation  analogous  to  (5).  An 
"adjoint"  scattering  configuration  (Figure  2)  may  thus  be  identified  from 
(16),  which  is  related  to  the  original  problem  by  interchanging  the  parameters 

“  -►  =  ss^r 

3.  «-»-  k.  •*-*•  -k  ,  and  U  •»->■  U  .  If  we  indicate  the  parameters  charac- 

terizing  the  original  field  E  explicitly  as  E(A,  kg,  k^,  Ip;*es,  3.;  U), 

then  these  interchanges  suggest  an  adjoint  field  solution  E  to  the  wave 
equation  for  the  configuration  in  Figure  2,  of  the  form  E  =  ..A*  -ki’  -ksi 
e.,  Ip*es;  UTr).24  By  deriving  the  adjoint  quantities  Nj,  N2,  and  D  for  E, 
analogous  to  (7)  through  (11),  one  observes  that  A  =  D/N2  =  D/Ni  since  D 
=  D  and  Ni  =  N2.  These  relationships  show  that  T  =  T  or  more  explicitly, 
£s-T(ks,  ki,  T-g$,  ai;  U)  =  g^T^,  -k$,  T?-g$;  UTr),  which  is  in 

the  nature  of  a  reciprocity  relation. 
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4.  VARIATIONAL  EXPRESSIONS  FOR  PERFECT  CONDUCTORS 

Analysis  of  scattering  from  systems  with  high  conductivities  is  often 
simplified  by  assunlng  the  scatterers,  e.g.,  rough  metallic  surfaces,  to 
be  perfect  conductors.2'25  In  this  section  we  shall  present  two  variational 
principles  for  perfect  conductors,  distinguished  by  the  type  of  dyadic 
Green  function  chosen  for  the  Integral  field  equations.  The  first  formula¬ 
tion  Is  specifically  applicable  to  rough  planar  scattering  surfaces.  The 
second  formulation,  advantageous  for  nonplanar  scatterers,  has  appeared 
elsewhere10  and  Is  Included  for  completeness.  As  for  the  previous  section, 
the  principles  derived  here  will  be  recast  Into  forms  suitable  for  stochas¬ 
tic  systems  In  section  5. 

4.1  Variational  Principle  for  Perfectly  Conducting,  Planar  Rough  Surfaces 

The  scattering  configuration  Is  Illustrated  In  Figure  3  in  which  a 
plane  wave  Is  Incident  upon  an  Infinite,  perfectly  conducting,  arbitrarily 
rough  planar  surface  So  which  we  will  restrict  to  lie  on  or  above  the  2=0 

plane.  In  the  linear,  homogeneous.  Isotropic  medium  (y,  e)  above  S0  the 

-► 

wave  equation  for  the  electric  field  E  is  the  homogeneous  vector  equation 

VxVxE(r)  -  k2E(r)  =  0,  (17) 

with  k2  =  u)2ye.  On  So  the  field  E  is  subject  to  the  boundary  condition 

4. 

nxE(r)  =  0,  reS0,  where  n  Is  the  outward  normal  from  S0.  On  the  surface 
at  Infinity  E  is  subject  to  the  radiation  condition  (3). 

To  formulate  the  Integral  equation,  we  choose  the  half-space  dyadic 
Green  function  of  the  first  kind  Gi,  which  satisfies  the  boundary  condition 
zxGi(r,r')  =  0  when  z  *  0  or  z'  =0.  This  function  has  the  form20 

G»(r,r')  =  {T  -  -p-  VV'j  [G(r.r')  -  G(r,rJ)]  +  2zz  G(r,ri),  (18a) 


where  V  and  V'  are  gradient  operators  over  r  and  r',  G(r,r')  =  exp(ikR)/4irR, 
and  the  Image  source  point  rj  is  related  to  r*  by  rj  =  r'  -  2z’z.  Equiva¬ 
lently,  G^r.r')  may  be  expressed  in  terms  of  the  free  space  dyadic  Green 
function  G  defined  In  section  3  as1* 

Gi ( r ,r ' )  =  G( r,r ' )  -  G(r , r • ) -R-  (18b) 
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/A 


In  which  R;  e  I  -  2zz  Is  the  reflection  operator  about  the  z  =  0  plane  (note 

rj  =  * r ' ) .  The  first  form  (18a)  exhibits  clearly  the  boundary  condition 

satisfied  by  6lt  and  the  second  form  (18b)  Illustrates  the  separate  field 
contributions  arising  from  a  source  and  Its  Image  In  the  presence  of  a 
planar  scattering  surface. 

By  proper  application  of  Green's  Theorem  to  (17)  In  the  volume  bounded 
below  by  So  and  above  by  the  surface  at  Infinity,  one  obtains  the  Integral 
equation 

E(r)  =?(?)+  Er(?)  -  lap  js  Gi (r,r*)-K(r')dS* ,  (19) 

■+  4 

where  K(r)  Is  the  surface  current  defined  by 

K(r)  =  nxH(r)  =  ^  nxVxE(r),  reSo, 

-4  -4 

-1  iki*r 

and  the  boundary  conditions  have  been  applied.  Here,  E  =  Ae..  e  is  the 

=  i (Rj’kj ) *r 

prescribed  Incident  plane  wave,  and  E  =  e  is  the  plane 

wave  which  would  be  specularly  reflected  from  a  perfectly  conducting  plane 
surface  at  z  =  0. 

Equation  (19)  is  valid  when  r  Is  not  on  the  surface  S0.  To  evaluate 
the  field  on  the  surface  [which  will  be  necessary  in  order  to  obtain  the 

Integral  equation  (26)  for  K],  care  must  be  exercised  in  view  of  the  singu¬ 
lar  behavior  of  Gi  at  R  =  0.  The  field  on  the  surface  is  defined  as  the 
limit  of  (19)  as  r  approaches  S0,  in  analogy  to  Reference  9,  Ch.  9.  This 

limit  can  be  verified  to  be  finite.  We  observe  that  the  choice  of  Gi 
results  In  integration  over  only  the  rough  region  of  So.  By  restricting 
the  surface  roughness  to  a  finite,  but  possibly  large,  portion  of  the 
Infinite  surface  S0,  the  asymptotic  form  of  (19)  becomes 

■+  ■*  +P  ■+  pi  kr 

E(r)  f^rrE1  +  Er  -  A  T  ,  (20) 

where 


T(ks,ki)  =  ^ftMk,.,^)  (21) 
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and 


Mk^)  • 

dS  I- 

_  -Ik  -r*  _  -1(Rj»k  )»r' 

I  e  s  -  R-  e  z  s 

4 

So 

. 

•K(r').  (22) 


The  amplitude  A  of  the  Incident  w§ve  may  be  eliminated  from  (21)  by 
first  Introducing.^  "adjoint  current",K(r)  which  we  assume  to  be  confined 
to  So,  such  that  K(r)  ■  0,  r^So  and  n*K(r)  ■  0,  reS0.  Premultiplying  (19) 
by  K(r),  Integrating  over  S0,  and  recognizing  the  boundary  condition  for 
E  on  So,  we  have 

0  -  |  dS  K(r)*E(r)  -  ANa(k$,kf)  -  1«w  D(ks,kf),  (23) 


where 


->  ^  -  1k.*r  _  1(Rc*k.)  r 

N«(ks,k1)  =  dS  K(r)«  21  e  1  -  R-21  e  z 

^So 


D^.k,)  s  f  dS  |  dS'  K(r)-G,(?,r’) •£(?•) .  (25) 

So  So 

In  accordance  with  the  discussion  following  the  Integral  equation  (19),  D 

in  (25)  Is  to  be  evaluated  by  first  performing  the  r’  Integration  with  r 

-*■  -*■ 

above  S0,  and  then  letting  r  approach  S0  to  perform  the  r  Integration. 
Combining  (23)  with  (21)  and  defining  the  projections  T  =  8  *f  and  Ni  = 
e.*Ni  yields  thg  familiar  form  (12),  which  Is  now  a  functional  of  the 
currents  K  and  K. 

To  Identify  the  conditions  for  variational  Invgrlance,  l.e.,  6T  =  0, 
we  perform  the„var1at1on  of  T  with  respect  to  K  and  K,  In  which, the  varia¬ 
tions  6K  and  6K  will  be  restricted  to  So,  with  n*6K  *  0  and  n*6K  =  0,  but 
are  otherwise  arbitrary.  The  condition  6T  *  0  yields  the  Euler  equations 


A  I-  81  e 


1k^*r 


1(Rj*k^)*r 


Rs*l<  e  =  luiy  dS '  I-*Gi(r,r')*I^,*K(r'),  reS0 
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and 


_  -Ik  *r'  _  -1(R;*k.)*r‘ 

I  e  s  -  R$  e  2  8 

•  I**£  * 

z 

r  s 

dS  I~.*Gi(r',r)*T-*K(r),  r'tSo 

(27) 


where  we  have  Introduced  the  projection  operator  I*  =  I  -  nn  onto  So  to 
emphasize  the  tangential  nature  of  the  currents  and  their  variations. 


[e.g. 


I-K(r) 


K(r)] 


Equation  (26)  Is  just  the  equation  for  the  current  K  on  S0,  as  we^ 

•f 

readily  deduce  from  (19),  while  (27)  Is  the  corresponding  equation  for  K. 
Conversely,  when  solutions  to  the  Euler  equations  (26)  and  (27)  exist,  then 

the  functional  T  Is  stationary  about  the  exact  currents  K  and  K.  Note 
that  the  left  side  of  (27)  may  be  Identified  as  the  tangential  projection 
onto  S0  of  the  Incident  and  specular  plane  waves  for  the  adjoint  problem 
shown  In  Figure  4.  As  In  the  case  of  scattering  from  conducting  dielectrics 

4  -►  *4  4 

In  section  3,  we  may  demonstrate  that  the  function  E  =  E(A,  -k^,  -ks>  e^, 

) »  where  ?(A,  ks,  k^,  I~*2s,  3.j)  Is  the  original  field,  satisfies  the 
adjoint  field  equation  and  Implies  ( Icay )~ 1  D/N2  *  (l^u)"1  D/Ni  *  A. 


4.2  Variational  Principle  for  Closed,  Perfectly  Conducting  Scatterers 

The  previous  "half-space"  variational  principle  Is  advantageous  for 
a  planar  rough  scattering  surface  because  the  Integrations  in  (12)  will  be 
non-zero  over  only  the  rough  portion  of  the  surface  S0.  However,  when  S0 
Is  closed,  a  scattering  expression  based  on  the  free  space  dyadic  Green 
function  (6)  Is  more  convenient. 10  The  variational  principle  for  the  con¬ 
figuration  In  Figure  5  has  the  standard  form  (12), 10  where  now 


Ni(ks,k1 ) 


j  -Ik  -r’  „ 

|  e  5  i$-I--K(r,)dS', 

Js, 


(28) 


i  -  1k,*r 

K(r)-ai  e  1  dS, 


and 


(29) 
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0(ks.k1) 


(30) 


|  dS  I  dS' 
Js0  'So 


K(r)‘G0(rtr,)«K(r,)1 


i 


The  Integrals  In  (30)  are  defined  by  the  same  limiting  procedure  used  In 

(25).  Straight-forward  variation  shows  the  Euler  equations  for  K  and  K 
to  have  the  form  (26)  and  (27),  but  without  the  explicit  appearance  of 
specular  terms. 

5.  STOCHASTIC  VARIATIONAL  FORMULATIONS 

The  variational  principles  of  the  previous  sections  will  now  be  recast 
Into  specific  forms  expressing  statistical  moments  of  an  arbitrary  component 

T  =  2$*T  of  the  scattering  amplitude  T.  Recall  that  each  of  the  variational 

expressions  has  the  general  form 


T 


1  NiN 
ti T  D 


(31) 


with  the  appropriate  Integrals  Ni,  N2,  and  D.  Following  the  procedure  of 
Hart  and  Farrell7  we  shall  use  the  nonstochastic  nature  of  the  amplitudes 
of  the  Incident  plane  waves  to  derive  exact  expressions  for  the  statistical 
moments  of  T  and  |T|2  In  terms  of  the  corresponding  moments  of  Ni,  N2,  and 
D.  More  Importantly,  the  expressions  will  be  proven  to  be  variatlonally 
invariant  as  a  direct  result  of  the  Invariance  of  the  expression  (31)  for 
T.7 


i 


We  begin  by  noting  that  for  each  scattering  configuration  equations 
of  the  form 


AT.^M,  (32) 

and 

A  N2  =  bD  (33) 


were  derived,  where  b  *  Iwy  for  perfectly  conducting  scatterers  and  b  =  1 

for  conducting  dielectrics.  In  addition,  our  choice  of  solution  E  to  the 
adjoint  wave  equation  In  the  previous  sections  led  to  ANi  *  bD.  If  the 
scatterers  fluctuate  randomly,  e.g..  In  shape,  position,  orientation,  or 
material  parameters,  the  expressions  (31),  (32),  and  (33)  remain  valid  and 
exact,  provided  that  we  use  the  exact  fields  or  currents  on  or  within  the 
scatterers.  Since  A  Is  the  nonstochastic  amplitude  of  the  Incident  plane 
wave  (and  since  b  Is  nonstochastic),  the  following  relations  may  be  obtained. 


(34) 


<Tn> / <N i n)  =  <Tn/Nxn>  =  Tn/N i n  -  bn/(47rA)n, 

(N2n>/<Dn)  =  <N2n/Dn>  =  N2n/Dn  =  bn/An,  (35) 

and 

<Nin>/<Dn>  *  <N2n/Dn>  =  N,n/Dn  =  bn/An,  (36) 

where  <  )  denotes  the  average  and  the  Integer  n  Is  arbitrary.  Combining 
(34)  and  (35)  results  in 


<T">  •  (^]n<N1n)<N2n>/<Dn>.  (37) 

This  is  an  exact  expression  for  <T  >  when  the  exact  fields  E  and  E,  or 
exact  currents  K  and  K,  are  Inserted  in  the  appropriate  integrals. 

We  will  show  that  if  6T  =  0,  then  6(Tn)  =  0,  that  is,  (37)  is  sta¬ 
tionary  about  the  exact  fields  or  currents.  First  define 


<T,n> 


((N1+6N1 )n)((N2+6N2)n) 
<(D+6D)n> 


(38) 


where,  for  example,  Ni  +  6Ni  is  the  result  of  approximating  the  exact  cur- 

rent  K  by  the  trial  function  K'  =  K  +  6K.  Thus,  the  trial  approximations 

for  K  and  K  (or  E  and  E)  introduce  the  errors  6Nj,  6N2,  and  6D,  which  are 
due  to  the  variations  of  these  currents  (or  fields)  from  their  exact  values. 
Then,  to  the  first  order  in  6Ni,  6N2,  and  6D, 


<5<Tn)/<Tn>  =  <T,n  -  Tn)/<Tn) 


=  n 


<(5Ni/Ni)Nin>  +  (( 6N2/N2)N2n)  _  < (6D/D)Dn> 


<Nsn) 


(N2°> 


<Dn> 


Applying  (35)  and  (36)  through  (39)  results  in 

6<Tn)/<Tn>  =  n<0n(6Ns/Nj  +  6N2/N2  -  60/D))/<Dn) 


(39) 


(40) 


For  variational  principles  with  the  form  (31),  the  quantity  In  parentheses 
in  (40)  vanishes,  and  we  have  proven  that 


6<Tn> 


(41) 
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when  Eq.  (36)  Is  used  to  express  (Tn>. 

To  calculate  statistical  moments  of  the  differential  scattering  cross 
section,  we  may  derive  a  general  variational  expression  for  (|T|2n).  Ap¬ 
plying  (32)  and  (33)  the  appropriate  expression  Is 


<|T|2"> 


2n 


<|Ni|2n><|N2|2n)/<|D|2n> 


Proceeding  as  before,  we  readily  show  that 


(42) 


<S<|T|2n>  *  0,  (43) 

that  Is,  first  order  variations  of  (42)  about  the  exact  fields  or  currents 
vanish. 

Under  the  usual  conditions  for  convergence,  note  that  the  Fourier 

•f 

transform  of  the  probability  density  function  f(T)  of  the  component  T  =  e  *T 
may  be  expressed  as  the  sum  of  statistical  moments,2*  s 

F(W)  =1  -^<Tn>,  (44) 

n=0  "• 

where  F  denotes  the  Fourier  transform  over  parameter  v.  A  similar  expres¬ 
sion  may  be  written  for  f(|Tl2)  where  [T|2  Is  simply  the  differential  scat¬ 
tering  cross  section  do/dfi.  Thus,  variational  approximations  for  f(T) 
and  f(|T|2)  may  be  obtained  from  the  variational  calculations  of  successive 
statistical  moments  of  T  and  |T|2. 

Variational  principles  of  the  form  (36)  or  (42)  are  readily  observed 
to  be  Inherently  more  tractable  than  corresponding  variational  expressions 
obtained  directly  from  (31),  l.e., 

<Tn>  =  (47)n<N,nN2n/Dn>  (45) 


and 


<min>  ■  (*j]!n<|N,l!"|NI|!rV|D|in>.  (46) 

We  see,  for  example,  that  (46)  Involves  calculating  the  average  of  the  com¬ 
plicated  quotient  NinN2n/Dn,  whereas  (36)  involves  the  quotient  of  the 
averages  <Nin>,  <N2n),  and  <Dn),  which  are  generally  simpler  to  evaluate. 
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6. 


SUMMARY 


The  configurations  for  which  stochastic  variational  expressions  of 
the  forms  (36)  and  (42)  have  been  derived  (Figures  1,  3,  and  5)  are  repre¬ 
sentative  of  a  wide  range  of  scattering  problems.  For  example,  precipita¬ 
tion  and  certain  aerosols  may  be  modelled  as  ensembles  of  conducting,  di¬ 
electric  scatterers,  to  which  (8),  (10),  and  (11)  are  Immediately  applicable, 
and  the  surface  of  the  ocean  may  be  modelled  as  a  rough  planar  conductor 
for  which  (22),  (24),  and  (25)  are  appropriate.  Stochastic  variational 
principles  show  promise  In  allowing  application  of  the  variational  method 
to  these  and  other  random  scattering  problems.  Future  work  will  Involve 
demonstration  of  their  tractablllty  by  explicit  calculation  of  the  statis¬ 
tics  of  electromagnetic  wave  scattering  from  models  of  random  volume  and 
rough  surface  phenomena. 
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APPENDIX 


FIGURES 


e  COMPONENT  OF  THE  SCATTERED  WAVE 
s 


Figure  I.  Geometry  for  Scattering  from  an  Arbitrary  Volume  Vo  with 
Permittivity  e(r)  and  Conductivity  o(r)  Suspended  in  Free  Space 
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Figure  2.  Counterpart  to  Figure  1  where  Roles  of  Polarization  and  Propa¬ 
gation  are  Interchanged  between  Incident  and  Scattered  Fields,  and 

e  and  a  of  the  Scattered  Volume  are  Transposed. 
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INCIDENT  WAVE 


es  COMPONENT  OF  THE  SCATTERED  WAVE 


V  k, 


Figure  3.  Scattering  of  a  Plane  Wave  by  a  Perfectly  Conducting  Planar 
Surface  with  a  Confined  Roughness  Region. 

(The  variational  principle  Is  derived  for  the  nonspecular 
portion  of  the  scattered  field.) 
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ei  COMPONENT  OF  THE  SCATTERED  WAVE 


nxE(r)  =  0  on  So 

Figure  4.  Adjoint  Counterpart  of  Figure  3  where  the  Roles  of  the  Incident 
and  Scattered  Wave  Parameters  are  Interchanged. 
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g  COMPONENT  OF  THE  SCATTERED  WAVE 
s 


INCIDENT  WAVE 


«i<  ' 


n*E(r)  =  0  on  Sc 


Figure  5.  A  Plane  Wave  Impinging  on  a  Closed,  Perfectly  Conducting  Sphere. 


F-j&'W  v  •'A  t  » •  -  'i  • 
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